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Abstract 



We present some analytic, non-perturbative results for the invariant density p{x) 
" ^ , for noisy one-dimensional maps at fully developed chaos. Under periodic boundary 

^1^ . conditions, the Fourier expansion method is used to show precisely how noise makes 

O ■ p{x) absolutely continuous and smoothens it out. Simple solvable models are used to 

(^ ' illustrate the explicit dependence of p{x) on the amplitude r] of the noise distribution, 

Q>^ \ all the way from the case of zero noise (r/ — > 0) to the completely noise-dominated 



limit (r/ = 1). 

PACS Nos. 05.45. +b, 05.40.+J 



1 INTRODUCTION 

One-dimensional (ID) maps exhibiting chaotic behaviour have been used as effective 
models to describe a wide variety of physical processes ranging from irregular behaviour in 
electronic circuits and chemical reactions to turbulence [1,2]. In particular, probabilistic 
or statistical approaches (involving, for instance, invariant measures on attractors) enable 
us to by-pass the limitations imposed by extreme sensitivity to initial conditions, and 
compute various quantities of interest in terms of statistical averages [3-5] . A considerable 
body of mathematically rigorous results on ID discrete-time dynamics is also available 
[6,7]. As random noise is inevitably present in physical systems, much effort has been put 
in to understand the effects of noise upon different aspects of chaotic dynamics, including 
invariant densities and related quantities [5, 8-17]. As much (but not all) of this work is 
based on numerical analysis or simulation, analytical results are of importance - especially 
so if they are non-perturbative in nature, i.e., if they are valid for arbitrarily large noise 
amplitudes. 

In this paper, we re-visit ID maps in the regime of fully-developed chaos, in the pres- 
ence of uncorrelated noise [18,19]. With the addition of noise, the evolution equation 
becomes a stochastic difference equation, thereby making the system effectively infinite- 
dimensional. Imposing periodic boundary conditions, the Fourier transform method 
[10,11] is used to establish some non-perturbative results for the invariant density p{x). 
In Sec. 2, we show how the addition of noise of arbitrarily small amplitude can make p{x) 
absolutely continuous in the interval. In Sec. 3, we consider a noise density that has a 
single non-trivial Fourier mode, and show that p{x) is also locked in at the same mode, 
but with a phase shift. In Sec. 4 we examine the dependence of p{x) on the amplitude 
of the noise density, varying the latter over the full range from zero (the noise-free limit) 
to unity (the completely noise-dominated limit). We conclude with a few remarks on the 
effect of noise on the Lyapunov exponent. 

2 CONTINUITY OF THE INVARIANT DENSITY 

We consider ID endomorphisms Xn+i = f{xn) , "^ = 0, 1, ■ ■ ■ , a;„ G [—1, 1] = I, exhibiting 
fully developed chaos with invariant density p^^\x). In the presence of additive noise, 

Xn+l = f\Xn) + in+1 (1) 

where (^„) = 0, {C,nC,n') = {^'^)Snn'- We use periodic boundary conditions, so that the 
normalized noise distribution g{C,) is also a periodic function with fundamental interval I. 



The invariant density p{x) for the noisy map (1) satisfies the perturbed Frobenius - Perron 
equation [5] 



p(x) = J^dyg{x- f{y))p{y) . (2) 



The Fourier expansions of p and g are 



1 OO -| oo 

P(^) = 2 ^ ^" exp{i7rnx) , g{^) = - J^ 9n exp{i7rn^) , (3) 



where 



Pn = L dxp{x) exp{-i7rnx) , gn = L d^ g{^) exp(-mO . (4) 

On substituting Eqs. (3) and Eqs. (4) and using the normahzation of p and g (^ 
Pq = 1, go = 1), Eq. (2) becomes equivalent to the following (infinite-dimensional) 
inhomogeneous matrix equation [10,11] for p„ , n 7^ : 



Pn 



g-n »-'nO ~r / ^ g-n ^nm Pm \^) 



where Y^' stands for a summation over all non-zero integers. The elements of the matrix 
S are given by 

Sum = - j dy exp[i-n(my-nf{y))]. (6) 

Iteration of Eq. (5) provides a fast-converging means of numerical solution, once the cjn 
are specified from the noise distribution concerned. The unperturbed (or noise-free) case 
is recovered on replacing g{^) by 5(^), i.e., by setting gn = 1 for all n. 

It turns out that exact solutions for the invariant density p{x) can be obtained from 
Eq. (5) in several important cases by exploiting the symmetries, if any, that the matrix S 
may happen to possess [17]. For instance, if Sno = for all n 7^ 0, both p'^'^^x) and p{x) 
reduce to the constant density 1/2. Again, suppose the noise density g{C,) is symmetric 
about its mean value zero (a reasonable assumption), so that ^„ = ^_„. Then, if Snm is 
even [respectively, odd] in the index m, while Sno is odd [respectively, even] in the index 
n, only the leading term in the iterative solution of Eq. (5) survives. This leads to the 
exact solution 

/'(^) = 2 ^ gnSnoexp{-i7rnx). (7) 

n=— 00 

Puting in the definitions of ^„ and Sno, Eq- (7) can be re-written in the form 

Pix) = - j^dyg{x- f{y)). (8) 



We shall use these results in Sec. 4 to study analytically the dependence of the invariant 
density on the amplitude of the noise distribution. 

Here, we wish to point out a simple way of understanding precisely how the addition 
of noise leads, in general, to a smoother invariant density [20,21]. Consider the noise-free 
case: The Fourier coefficients of p^^' (x) are given by 

m m ml 

The asymptotic (large n) behaviour of p^°) is thus controlled by Snm- If Pn^ ~ n^^, then 
p^^\x) has finite discontinuities in I - including, possibly, the end points ±1, as we have 
used periodic boundary conditions. Now consider what happens when noise is added. If 
g{^) is continuous, then Qn ~ n^"^. Even ii g{^) has finite jumps in I (including, possibly, at 
the end points ±1), its aymptotic behaviour is at least 0{n^^) [22]. It follows at once from 
Eq. (5) that the asymptotic behaviour of p„ is improved to p„ = 0{n~^''^) where a > 0. 
Consequently, J2 Pn exp{i7Tnx) converges absolutely in I, and p{x) becomes continuous 
everywhere, including the end points. (That is, p(— 1) = p(l)-) Explicit illustrations will 
be given subsequently. 

3 EXACT SOLUTION FOR 'SINGLE - MODE' 
NOISE DENSITY 

We show now that an interesting form of "mode-locking" occurs if the noise density has 
a single wavelength, i.e., Fourier component. The corresponding noise density is given by 
the one-parameter family of functions 

9(0 = ^[1 + (-ir'cosTrr^], (10) 

where r is an integer. Correspondingly, 

1 

gn = 5nfl + -{-iy''^{5n,r + 5n-r)- (H) 

Hence the only nonvanishing Fourier coefficients p„ {n ^ 0) are pr and p_,. = p*. Using 
the relation S'_„,_m = •S'* „, we find from Eq. (5) the solution 

0( 1 V' — ^ Q c C* _i_ C* C* 

~ _ ^l -"-J '-'rO 'JrO 'Jrr "•" ^rO '-'-rr t T 0\ 



Writing pr = \pr\ exp(i0), this leads to the result 



P(^) = 9 + \Pr\ cos{-Krx + (j)) . (13) 



The invariant density is therefore 'locked in' at the same wavelength as the noise density, 
with a phase shift. The amplitude and phase shift depend on the map / through the 
matrix elements of S that appear in Eq. (12). It is evident that the situation corresponds 
to completely noise-dominated dynamics. This is brought out explicitly in the example 
considered in the next section. 

4 DEPENDENCE OF p{x) ON THE NOISE AMPLI- 
TUDE 

We want to study the explicit dependence of the invariant density on the amplitude of 
the noise density in a non-perturbative manner. For this purpose, we consider illustrative 
cases in which an exact solution for p{x) is possible owing to symmetries present in Snm, 
as indicated in Sec. 2. The noise density is taken to have a compact support [—Ti,ri] (i.e., 
an amplitude r]). A convenient form for the normalized density g{^) that enables us to 
scan the entire range of rj from to 1 is 

f (2ri)-^(l + cos^), 1^1 < r7 
I 0, V<\^\<1 

where < r^ < 1. Figure 1 depicts g{^) for rj = 0.2. In the limit rj ^ 0, we have 
g{C,) -^ ^(0) '^^ noise- free dynamics. When t] ^ 1, we have the fully noise dominated 
case considered in the preceding section. The Fourier coefficients of g{^) are given by 

smimrr]) \ -n? rf 1 
Qn = 1 ^^ • 15 

To bring out all the effects of noise on p{x) which we want to demonstrate, let us first 
consider the square-root cusp map (see Fig. 2) 



Xn+l 



f{Xn) = l - 2\Xn\^/^ , X„el. (16) 



It is well known that this map provides a prototypical model of intermittent chaos, arising 
from the marginal stability of the fixed point at the left boundary, x = —1. However, 
unlike many other models of intermittency that share the latter feature, the map (16) has 
the non-singular invariant density [23] 

p^'\x) = i^ . (17) 

As f{x) is symmetric, we have Snm = Sn-m- Further, 

1 /• . .. . „ (-1)" 



1 f (~lj 
Sno = - dy exp{2Tiiny/y) = — — - = -S. 

2 Jl nirn) 



n,0 • 



Hence, as pointed out in Sec. 2, the exact solution for p{x) is given by Eq. (7). We get 

1 °° (—1)" . 

P{^) = o + H gn sin^nnx) . (19) 

2 „=i nn 

Putting in the expression for cjn from Eq. (15) and carrying out the summations involved, 
we finally obtain the closed-form expression 

p{x) = { i^, -i^r^<x<l-r] (20) 



I- ±(l-r])il-x) + ^sin^(^, l-r]<x<l. 
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Figure 3 shows p{x) in the cases rj = 0.1, 0.2 and 0.7 respectively. (Throughout this 
paper, relatively large values of r] have been chosen for clarity of illustration, and also to 
emphasize the non-perturbative nature of the results.) The following points are notewor- 
thy: 

(i) Both p^'^\x) and p{x) are antisymmetric about the mean value 1/2. However, p^'^\—l) 
= 1 while p*-°-'(l) = 0. But p{x) is continuous everywhere, including the end points 
±1 (p(±l) = 1/2), in accordance with the general result established in Sec. 2. This 
feature persists for arbitrarily small values of rj. It is known that it is rather difficult in 
numerical simulations to obtain the exact result p^^\x) = (1 — x)/2: a 'boundary layer' 
persists, especially near a; = — 1, in which the invariant density builds up to a value close 
to unity, instead of straightaway starting at that value for x = —1 and then decreasing 
linearly as x increases. Our result on the continuity of p{x) accounts for this phenomenon: 
noise, albeit in the form of round-off errors, is inevitably present in simulations, (ii) In 
the opposite, noise-dominated limit r^ ^ 1, we find 

1 sin (irx) 

again in accord with the results of Sec. 3, in the special case r = 1, \pr\ = l/(27r) , = 
7r/2. 

Next, let us consider a simple case in which p^'^\x) has a finite jump in the interior of 
the interval, in order to see how the noise smoothens it out. A convenient illustration is 
provided by the piecewise linear map (see Fig. 4) 

J 1 ~ 5|x| , |a;| < I 

The unperturbed invariant density is the piecewise constant function 

rm. s f t, -1< a; <0 
i |, 0<x <1 
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with a jump at a; = 0. In this case, too, Snm = Sn,-m, and 

5nO = -S-nO = ^[1 - (-1)1 . (24) 

When noise (distributed as in Eq. (14) is added, p{x) is found to be the following piecewise 
analytic function, for r^ < 1/2: 

—l<x<—l+ri 

— 1+1]<X< —1] 

Pi^) = S i - ^ - T^ sin if , \x\<V (25) 

7] < X < 1 ~ f] 
1 ~ 1] < X < 1 . 

p{x) remains piecewise analytic for rj > 1/2, as well. As before, we illustrate in Fig. 5 
the trends in the variation of p{x) using large values of rj : 0.2, 0.5 and 0.7 respectively. 
As expected, the discontinuities that p^^\x) has at a; = and ±1 disappear in p{x), for 
arbitrarily small rj. As rj ^ 1, p{x) tends to 1/2 — (S/Stt) sin(7ra;), corroborating our 
result obtained for 'single - mode' noise. 

5 CONCLUDING REMARKS 

We have shown that additive noise significantly alters the invariant density for fully chaotic 
ID maps: arbitrarily small amounts of noise smoothen out the density and remove any 
discontinuities it may have in the absence of noise, under periodic boundary conditions. 
One may ask whether the Lyapunov exponent A*-^-' is also altered by the addition of noise 
[24]. While this is so in general, there is no change in A'-'^-' in the examples considered in 
Sec. 4, for the following reason: if the map f{x) is an even function of x, so is In |/'(a;)|. 
If, further, both p'^'^^x) — 1/2 and p{x) — 1/2 are antisymmetric in x, it follows that the 
corresponding Lyapunov exponents are equal to each other, and are given by 

A(°) = a = / dx\n\f'{x)\. (26) 

Jo 

The general conclusion is as follows: if f{x) is an even function in I and so is the 
noise density, then A remains equal to A'-'^^ provided S'„o is odd in the index n, i.e., 
/j dx cos (nnf^x)) vanishes for every integer value of n. 
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Figure Captions 

1. Plot of^(0 vs^hiri = 0.2 

2. The square-root cusp map 

3. Plot of p{x) vs X for the noisy square-root cusp map corresponding to rj = 0.1, 0.2 
and 0.7 respectively 

4. The piecewise linear map 

5. Plot of p{x) vs X for the piecewise linear map corresponding to ?7 = 0.2, 0.5 and 
0.7 respectively 
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